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ABSTRACT 


Computer models have been developed for analyzing the state-of-s tress 
beneath traction-drive type of contacts. The analyses involve computing 
stresses and stress-reversals on various planes for points beneath the sur- 
face. The effect of tangential and axial friction under gross slip conditions 
can be evaluated with the models. Evaluations performed on an RC (rolling 
contact) tester configuration indicate that the classical fatigue stresses are 
not altered by friction-forces typical of lubricated contact. Higher values 
of friction (f > 0.25) can result in surface shear reversal that exceeds the 
.'tresses at the depth of maximum shear reversal under rolling contact. 



STRESS EVALUATIONS UNDER 
ROLLING/ SLIDING CONTACTS 

by 

J. W. Kannel and J. L. Tevaarwerk (Consultant) 


INTRODUCTION 


Rolling elements such as roller bearings and traction drive trans- 
missions are subject to failure from rolling contact fatigue. Such fatigue 
failures cause a serious restriction on the operating life and reliability of 
such devices. In the case of the traction drive transmissions, rolling con- 
tact fatigue failures could conceivably restrict the use of such transmissions 
in motor vehicles, despite the increased efficiency such transmissions could 
afford. Although a number of factors contribute to fatigue failure, the ef- 
fect of traction and slip on the state of stress between rolling elements is 
not understood. A better understanding of the role of traction on fatigue 
could greatly assist in the evaluation of materials for and operational limits 
of traction components. 

The fundamental theory for rolling contact fatigue is the one pub- 
lished by Lundberg and Palmgren^)* in 1947 and summarized by Coy, et 
al.(2) # This theory is based on the assumption that failure is in the form 
of shallow pitting and is related to subsurface shearing stresses within the 
rolling elements. In essence, any irregularity in the material beneath the 
surface can manifest itself as a failure Initiation point. The significant 
stress field has been hypothesized as being bounded by the surface, the depth 
of the maximum and reversing shear stress, and the width of the rolling track. 
The magnitude and depth of the reversing shear stress have been found to be 
very significant with regards to fatigue failure. 

The reversing stress field has been well established for pure rolling 
contacts(3,4) , However, when surface tractions are imposed, this stress 
field becomes very complicated. The classic study for line contact is by 
Smith and Liu^). Smith and Liu's research Indicated that friction is a key 


^References are listed at the end of the text 
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factor with regards to reversing stresses. For example, at some level of fric- 
tion, the reversing stress was found to be at the surface, which would greatly 
reduce fatigue life. Hamilton and Goodman^) developed a theory for a cir- 
cular sliding contact. The results of their theory are in the form of useful 
graphs for predicting friction effects for this specific type of contact. 

Other researchers, such as Kuznetsov^), have studied the influence of fric- 
tion on contact stress, but not on the subsurface stress field under slip 
conditions. 

One complication to the reversing stress field is that the stress 
reversals occur over the width of contact and depend on orientation as well as 
depth. See, for example. Figure 1. The maximum stress occurs at the center 
of contact and the depth and magnitude are straightforward to compute. Com- 
puting the reversing stress becomes more complicated because the reversals 
occur over a finite region. Consider, then, the complication when the magni- 
tude of t + or x — depend on the orientation out of the plane of the paper. 

The goal of the project has been to develop computer models to deter- 
mine the magnitude of reversing shear stresses beneath (and very near) the 
surface of rolling/sliding contacts. Sliding in both the tangential direction 
and the axial direction are considered. Considerations are given for line 
contacts as well as crowned rollers. In order to insure accuracy in modeling, 
two independent models were developed, one by Tevaarwerk (T-model) and one at 
Battelle (B-model)*. This approach allowed for both researchers to be heavily 
Involved with the problem in modeling and created excellent dialogue, both for 
the mathematical intricacies, as well as the ramification of the predictions. 

A comparison between the two models is given in Table 1, which demonstrates 
the consistency of the two approaches. The approach for the development of 
the models was as follows: 

• Define the stress tensor for any point beneath the surface 

• Compute principal stresses and their direction cosines for each 
point 

• Determine the direction cosines (relative to the x, y, z axes) to 
the octahedral plane and the plane of maximum shearing stress 


*As discussed In Appendix A, the primary difference in the two models is in 
the Integration scheme used to determine the stress tensor. 
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FIGURE I. ILLUSTRATION OF SHEAR STRESS BENEATH THE SURFACE OF CONTACT 
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TABLE 1. COMPARISON OF TF.VAARWERK MODEL WITH BATTELLE MODEL 
(K - . 5 , v - . 285 , f A = 0 , z/b « . 3 ) 





°x 

a 
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°z 

T xy 
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I'* 
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f T 
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.228 
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-.256 

-.255 
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• Scan the planes of octahedral and maximum shearing stress in the 
direction of rolling and compute stress reversals. 

Stress reversals can be computed by the above procedure for condi- 
tions of zero friction, axial friction, and tangential friction for line 
contacts and contacts corresponding to a rolling contact fatigue tester. The 
basic assumptions for the analyses are: 

1. The contact pressure (normal stress) are predictable by Hertz 
theory for line or crowned contacts. 

2. The contact tractions are a result of a constant traction 
coefficient times the pressure. This tacitly assumes a gross 
sliding situation. 

3. The materials are all elastic and homogeneous. 

4. The bodies are isothermal. That is, no thermal stresses are 
considered. 


MODEL DEVELOPMENT 


The Stress Tensor 


The objective of the analyses is to evaluate subsurface stresses in 
rolling/ sliding contacts of the type seen in traction drives. Consider, for 
example, the system shown In Figure 2, where two cylinders are loaded together 
under rolling or possibly rolling/sliding contact. The assumed Hertzian con- 
tact pressure distribution produces a state of stress throughout the cyl- 
inders. Likewise, the frictional forces associated with sliding produce an 
additional state of stress such that a typical point beneath the surface 
is undergoing a complicated and changing stress state. For crowned cylinders, 
this stress state varies axially (in the y direction) as a result of axially 
varying pressures as well as in the x direction, as shown in Figure 3. 
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In tensor notation, the stress at any point can be written as 

°ij ° r , a) 

where °ij ” a lj/PH (i, j "* 1, 2, 3) , 

where is the maximum Hertz contact pressure. In more conventional 
terminology 

°x ■ °11 * °y * °22 . °z * °33 • ^ 

T xy “ a 12 * T xz “ a 13 an< * T yz “ a 23 • 

It should be noted that the stress tensor is symmetric such that o^j • 

To determine o^j requires that expressions be developed for the 
effect of point loads on stresses and that the point loads be Integrated over 
the conditions of interest. These equations are discussed in Appendix A and 
Equation A-27 shows the general relationship for o’jj of Interest here. 

Principal Stresses 

Prom the preceding discussions. It can be noted that the state of 
stress at any point has nix components*. The components discussed in the pre- 
ceding section describe the stress when the axis Is oriented In the x, y, * 
system. However, the stress can look considerably different at different ori- 
entations. Since the objective of the study is to evaluate "worse case" shear 
stresses, it is necessary to evaluate the effect of axis rotation (see Figure 
4) on stress. The first step In this phase of the study is to determine the 
parameters that are not a function of rotation; these are known as the Invari- 
ants of the stress tensor^®). These Invariants are: 


*Actually nine components are required, but some are equal because of 
symmetry. 
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l r ■ °11 + °22 + °33 > ("0 

II r * °22 + a ll a 33 “* 0 22 a 33 ~ a 12^ 0 1 3^ ~ a 23^ * and 

1 1 1 r * a ll a 22 a 33 + 2 a\2 a 23 °13 ~ a ll a 23^ “ a 2 2 a l 3 ^ ~ a 33 a 12^ * 

i 

Z 



i 

FIGURE 4. ILLUSTRATION OF AXIS ROTATION 

With Che knowledge of the inverients, It le possible to determine the 
plane of the principal stress (i.e., the rotated axis on which no shear forces 
exist) and the stress aagnltudes. Tht equations for the principal stresses 
can be written 


S3 ♦ I r s2 ♦ H r S ♦ III r - 0 


( 6 ) 
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Equation 6 can be readily solved by cubic equation solutione and will produce 
three real roots for principal stresses. These roots will be known as , 

S2, and S3 with Sj being the largest and S3 the smallest. 

Angles To The Principal Stresses 

In tensor computations, the angles are given in terms of direction 
cosines to the various axes* If we define a nx , a n y, and a n2 as the 
direction cosines from a given principal stress (&ay S\) to the x, y, and z 
axis, then it can be shown that the following equations must be 

satisfied^) : 


a nxl^ + a nyl^ + a nxl^ m l » (7) 

(°11 ~ Si) a nx i + 0i2 a nyl + °13 a nzl * ^ 

°12 a nxl + < a 22 " s l) a nyl + a 23 a nzl " 0 » and 

a 13 a nxl + a 23 a nyl u <°33 * s l) a nzl * 0 • (10 > 


By solving Equation 7 and any two of Equations 8 through 10 simul- 
taneously, a nx i, a n yi and a nz i can be determined. In a ^ike manner, the 
direction cosines for locating the S 2 ( a ux2» a nz2» a nz2 ) and the S 3 
( a nx3 , a n y 3> aj,^) stress can be determined. 

Other Planes of Interest 


The solution of equation 7 through 10 permit the determination of the 
location of principal stresses. However, the objective of the project is to 
evaluate shear stresses and by definition there are no shear stresses on these 
planes. Therefore, other planes must be considered. 
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Plane of Maximum Shear Stress 


The maximum shear stress for an element lies on a plane bisecting the 
Si - S3 plane. The direction cosines for the normal to this plane rela- 
tive to the plane of the principal stress axes can be written 


b n2 = 
b n3 “ 

To relate these cosines back 
transformation of the form 

c nx = b nl a nxl 
c ny 3 b nl a nyl 
c nz 3 b nl a nzl 


1/S2 , on 
0 . , ( 12 ) 
1//T . (13) 
to the x, y, z axis requires an angle 

+ b n2 a nx2 + b n3 a nx3 » (14) 
+ b n 2 ^ny 2 + b n3 a ny3 » (15) 
+ b n2 a nz2 + b n3 a nz3 » ^6) 


where , c n 2 and c n 3 are the direction cosines from the normal (to the 
shear plane) to the x, y, z axes. 


Plane of Octahedral Stress 


The location of the normal to the octahedral stress can be written 


bnl " 1/^3 , 

(17) 

b n2 * 1/VT , 

(18) 

b n3 - 1//3 . 

(19) 


That is, the normal is situated at equal angles from the principal axis. 
Likewise, Equations 17 through 19 can be used in conjunction with Equations 14 
through 16 to relate the normal to the x, y, z axes. 



Seeking The Maximum and Minimum Shear Stresses 


Using the angle transformation equations discussed in the previous 
section, it is possible to locate the normals to various critical shear stress 
planes. The final step in this process is to locate maximum and minimum shear 
stress on those particular planes for various x positions. Let us define the 
direction cosines to a given location on a shear plane as a sx , a 8 y and 
a 8Z . It can be shown that^®) 


a sx c nx + a sy c ny + a sz c nz * ® 

and 

( 20 ) 

a 2 + o 2 + a 2 » i 

a sx + a sy a 8z 1 


( 21 ) 


If, then, we choose a particular angle, say a gx , we can compute 
a S y and a 8z by Equations 20 and 21. 

With a knowledge of the angles between the shear vector on a plane 
and the x, y, z axes and a knowledge of the stress tensor as oriented to the 
x, y, z axes, the shear stress can be computed. The equation for the stress 
is given by(®) 


T " c nx a sx °11 + c ny a sy °22 + c nz a sz 0 » 

+ c nx a sy °12 + c ny a sx °12 + c nx a sz °13 » ^2) 

+ Cnz a sx °13 + c ny a sz °23 + c n * a sy °23 • 

In the computer model, values of a 8X were assumed in the range 

-1 to +1. For each value of a sx , the shear stress was computed at various x 
positions and both the maximum and minimum values were stored. In the compu- 
tation, the maximum shear stress and the maximum reversal t ux - T a i n were 
sought. 
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Construction Of Model 


The computer model is constructed as follows: 

1. Input values of X (b/a), v (Poissons ratio), z/b (depth) f-p 
(tangential friction coefficient), f^ (axial friction 
coefficient), and case number. 

2. Compute the stress tensor for approximately eleven (11) x 
locations (STRSSC). For K = 0 (line contact), the Smith-Liu^^ 
solutions are used. 

3. Compute the principal stresses at each x location (PRINC). 

4. Compute the angles to the principal stresses at each location 
(ANGLP) . 

5. Compute the angle to the normals of the maximum shear plane, the 
octahedral shear plane, and the orthogonal plane (ANGTRAN). 

6. Compute the shear stresses and shear reversals on the shear 
planes (RVSTRS), 

All computations are in the form <?ij/PH and are applicable for 
any (elastic) contact pressure level. A source listing for the model is given 
in Appendix B and a typical printout is given in Table 2. 

In this table 
K * b/a (aspect ratio) 

NU * (Poissons ratio) 

Z * z/b (depth into surface) 

Fj » f-j. (tangential friction coefficient) 

^A * ^A (axial friction coefficient) 

Sjq(, Syy, Sg£ are the normal stresses (o x , 0y, o 2 ) 

S XY* S YZ» S XZ are the shear stress (r xy , t yz , t X z) 

Si, S 2 » S 3 are the principal stresses 
X - pos » x/b is the x position being evaluated 

TAU-MAX » maximum shearing stress as computed by scanning different 
planes 

AN (1), AN (2), AN (3) are the direction cosines to the normals on 
the planes containing TAU-MAX and, AS (1), AS (2), AS (3) are the 
direction cosines to TAU-MAX on the shear plane. 


TABLE 2. TYPICAL PRINTOUT FOR BATTELLE STRESS MODEL 


O 

o 


■ 

ft 


® 

z 

3 


o 

o 


o 

o 


■ 

3 


o 

o 

♦ 

UJ 

o 

o 

o 

o 


X 




rH 

O 

o 

o 

o 

O 

o 

O 

o 

rH 




o 

o 

o 

o 

o 

© 

© 

o 

o 

o 




4 

4 

♦ 

4 

♦ 

♦ 

♦ 

4 

♦ 

♦> 




UI 

ui 

UJ 

Ui 

UI 

ui 

ui 

UJ 

Ui 

UI 




o 

O 

o 

o 

o 

o 

o 

o 

o 

o 




o 

O 

o 

© 

o 

o 

o 

o 

o 

o 



A 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 



a 

o 

o 

o 

o 

o 

o 

© 

o 

o 

o 



QL 

rH 

© 

® 

A 

pp> 

PH 

A 

® 

o 

pH 



X 

r 

1 

1 

•* 

1 o 









o 

© 

o 

o 

© ^H 

o 

o 

o 

o 

© 




o 

o 

o 

o 

o o 

o 

o 

o 

o 

© 




♦ 

♦ 

♦ 


♦ ♦ 

♦ 

♦ 

♦ 

♦ 

♦ 




Ui 

Ui 

Ui 

Ui 

UJ Ui 

Ui 

UJ 

Ui 

Ui 

UJ 



ih 

r 4 

lA 

m 

A 

PM ® 

rH 

A 


o 

pH 



to 

PH 

4 > 

PH 

A 

A o 

PM 

o 

o 

Hi 

® 




PH 

Hi 

lA 

® 

H- Hi 

PH 

o 

r- 

pH 

o 




Hi 

PH 

® 

A 

A O 

O 

rH 

o 

P- 

o 




Hi 

4 

IO 

® 

O pH 

o 

o 

4 ) 

A 

IH 




• 

• 

• 

• 

• • 

• 

• 

• 

• 

» 




1 

1 

1 

1 

1 1 

1 

1 

1 

1 

1 




O 

O 

o 

O 

o o 

© 

o 

O 

O 

O 




o 

o 

a 

o 

o o 

o 

o 

o 

o 

© 




♦ 

♦ 

♦ 

♦ 

♦ ♦ 

4 

4 

♦ 

♦ 

♦ 



HI 

Ui 

Ui 

Ui 

UI 

Ui Ui 

UI 

u. 

UJ 

UJ 

Ui 



a 

Hi 

4 

*H 

PH 

pw g) 

4 > 

+ 

pH 

pH 

Hi 




<o 


pH 

® 

<r ® 

PH 

A 

4 

O 

fH 




HJ 

HJ 

A 

pH 

(H A 

rH 

A 

4 

rH 

43 




r- 

43 

4 

H 

rw m 

Hi 

A 

O 

pH 

fH 




rH 

HJ 

HI 

A 

A 4 ) 

4 ) 

A 

A 

PH 

Hi 




• 

• 

• 

• 

• • 

• 

• 

• 

• 

• 




1 

1 

1 

1 

1 1 

1 

1 

1 

I 

1 




pH 

H 

O 

© 

O O 

O 

O 

O 

o 

rH 




o 

o 

© 

O 

o © 

o 

© 

o 

o 

© 



rH 

1 

1 

♦ 

♦ 

♦ ♦ 

♦ 

♦ 

♦ 

♦ 

1 



A 

Ui 

Ui 

UI 

Ui 

Ui Ui 

UI 

Ui 

UJ 

UJ 

UI 




pH 

4 “ 

© 

© 

Hi PH 

o 

Hi 

4 ) 

A 





4 ) 

43 

pH 

A 

H* pH 

® 

K 

O 

A 

H- 




m 

O' 

A 

f- 

O PH 

4 

4 > 

rl 

>o 

® 




A 

Hi 

© 

A 

4 ) pH 

O 

4 

Hi 

® 

A 




Hi 

H* 

pH 

4 * 

A 43 

A 

A 

PH 

rH 

<n 




• 

• 

• 

• 

a • 

• 

• 

• 

• 

• 





1 

1 

1 

1 1 

1 

I 

1 

1 

1 




O 

O 

rH 

pH 

pH pH 

O 

O 

O 

o 

o 




O 

o 

O 

o 

o o 

o 

O 

o 

o 

o 





♦ 

1 

1 

1 1 

♦ 

♦ 

♦ 

♦ 

♦ 



PU 

UJ 

Ui 

UI 

Ui 

UJ Ui 

Ui 

Ui 

UJ 

Ui 

UJ 

>- 


X 

m 

fH 

4 “ 

o 

A A 

® 

m 

A 

o 

4 ) 

ft 


a 

® 

Hi 

HJ 

H- 

O' O 

o 

4 ) 

PM 

o 

O 

4 



in 

o 

4 

H- 

Hi O 

® 

A 

® 


4 

E 



rH 

rH 

H. 

O 

A pH 

o 

Hi 

43 

® 

H* 

Z 

• 


rH 

fJ 

H- 

pH 

4 - ® 

pH 

pH 

pH 

pH 

pH 

3 

o 


• 

• 

• 

• 

• • 

• 

• 

« 

• 

• 

A 






1 

1 1 

1 

1 

1 

I 

1 


a 












A 

x 


® 

® 

O 

® 

« ® 

® 

® 

OP 

® 

® 

A 

u_ 


o 

o 

o 

o 

o o 

o 

o 

o 

o 

o 

UJ 

i 


1 

1 

I 

1 

1 1 

1 

1 

1 

1 

1 

ft 

o 


UJ 

Ui 

UJ 

UJ 

Ui Ui 

Ui 

UJ 

Ui 

Ui 

UI 


o 

Hi 

A 

pn 

A 

4 

O 4 ) 

43 

4 

A 

PH 

A 

A 

♦ 

> 

A 

H 

© 

4 

4 * O' 

4 

4 

O 

pH 

A 


UJ i 

A 

PH 

Hi 

HI 

Hi 

4 - pH 

4 

PSJ 

*M 

Hi 

PH 

ft 

O 


A 

IA 

A 

A 

A A 

A 

A 

A 

A 

A 

4 

o 


ph 

PH 

fPS 

m 

A PH 

PH 

PH 

PH 

PH 

PH 

Ui 

o 


• 

• 

• 

• 

• • 

• 

• 

• 

• 

• 

X 

o 












A 

*4 













• 1 


<r 

« 

® 

® 

® « 

® 

eo 

A 

CD 

® 

UJ 



o 

O 

© 

© 

O © 

O 

O 

O 

o 

O 

3 



1 

• 

1 

1 

1 1 

1 

I 

1 

• 

1 

© 

a 

>- 

UJ 


J 

Ui 

UI UJ 

Ui 

Ui 

UJ 

Ui 

UJ 

M 

t- 

X 

A 


A 

4 

4 > 4 ) 


4 

A 

PH 

A 

K- 

IL 

A 

A 

w* 

© 

4 

4 - O' 

4 

4 

O' 

rH 

A 

4 



<H 

PM 

HJ 

Hi 

4 > p* 

4 

Hi 

PSi 

Hi 

PH 

U. 

O 


IA 

A 

A 

A 

A A 

A 

A 

A 

A 

A 


O 


PH 

ph 

PH 

A 

m PH 

PH 

PH 

PH 

PH 

PH 


♦ 


• 

• 

• 

• 

• • 

• 

• 

m 

• 

• 


Ui 













o 













o 


o o o 

OOHO 

o 

o o 


o 


o o o 

O O O O 

o © 

o o 


o 

(HI 

♦ 

♦ 

♦ 

♦ 

♦ ♦ 

4 

♦ 

♦ 1 

1 ♦ 

♦ 


rH 

rsi 

u> 

UJ 

UI 

Ui 

Ui Ui 

Ui 

UJ 

ui 

1 UJ 

Ui 


• 

A 

4 > Hi 

A 

<o 

o (> s 

-1 4 K 

43 Hi 




<0 O 

PH 

1 (Si PSi Hi A 


H- O 




® o © 

i ® Hi O ® 

O 

PH* 

O H» 


a 


PH © 


1 A A O A 

4 ) O 

pH 4 ) 


tej 


«-« 

in 

A 


O m» 

* 

® 



rH 




13 


a 

ft 


x 

a 


in O 

ft o 
I ♦ 
x ui 
PSJ 
O *H 
O lA 

♦ o 


cr» | 

ro 

h- 

43 a 
• <■— 
I <*> 


~ O 

A o 

W + 

Z Ui 

4 o 
o 

O M 

o o 
♦ f- 

UJ • 
H | 


o a 

r- 

• hi 

t/> 


«■* O 
rvj O 
— 4 
Z Ui 

4 O 

o 
o o 

0 o 

+ rH 

Ui • 
PH t 
® 

A 

rH • 
Hi -*> 

• H 

1 ^ 

A 

4 

« 

m O 

*-• o 
~ ♦ 
Z Ui 
<4 Hi 
O 
O rH 
O Hi 
♦ PH 
Ui • 
Hi 


O « 

Hi O 


• UI 

X > 


3 3 

4 4 


A 

A 


X 

a 


o 

o 


o 

o 

o 

o 


UJ O 

ft o 

I ♦ 

X UJ 

o o 

O HI 
♦ ® 
uj A 
4 > • 

H* | 

4 

o 

r- • 


a 

a x 

o 
a o 
— + 
Z U 
4 Hi 
4 
O O 

o 4 
♦ <o 

Ui • 
A 


4 

* O 

t o 
> ♦ 


Hi o 
Z ' 


[ o 
o 
o o 
o o 

♦ tH 
Ui • 
A | 


— o 

rH O 
♦ 


o 4> 
O H“ 
♦ Hi 


O 

o a 

HJ O 


o • 
o o 
o 
o 


A 

a 4 


a 

• HI 


o o 
♦ m 


© 

© 


«j a 

** © 


33 

4 4 


ELAM Of HEPEBS1NG SHEAR STRESS 






14 


TAU-REVERSING is the maximum reversing shear stress on the planes 
specified . 


* 


ANALYSIS FOR RC FATIGUE TESTER 


One objective of the project has been to evaluate the effect of mag- 
nit^de and direction of slip on the stresses in a rolling contact (RC) fatigue 
tester* The RC tester consists of a rotating cylindrical specimen loaded 
between two large diameter crowned rollers* In the test, the specimen is 
driven by an electric motor; the crowned rollers are only driven by the speci- 
men. The rollers are loaded to a level to produce fatigue in the specimen in 
a reasonable time period. Normally, the RC tester is used to evaluate rolling 
contact fatigue and has been used in bearing-material evaluations. The tester 
is also being considered as a candidate device for evaluating the effect of 
friction on fatigue. 

Two approaches could be used to induce friction between the 
cylindrical specimen and the crowned rollers. One technique would be to 
induce a drag between the two elements by driving or braking the rollers. 

This approach could be quite complicated for simple fatigue testing and an 
alternate could be more desirable. This alternate approach could consist 
simply of skewing the specimen relative to the rollers. However, skewing 
produces an axial traction, whereas braking the rollers produces a tangential 
traction. The stress model can be used in determining the stress fields 
associated with these two types of fatigue test concepts. 

Typical dimensions for an RC rig are 
D s - 9.5 mm (0.375 inch) 

D r * 190 mm (7.5 inches) 

R c - 6.4 mm (0.25 inch). 

D 8 is the specimen diameter, D r is the roller diameter, and R c is the 
crown radius. For these dimensions, a value of K * 0.8 has been computed. 

The upper limit traction coefficient should be nominally 0.1 for lubricated 
tractions. 

Figure 5 shows the effect of traction on the reversing shear stresses 
in the RC tester. As can be observed, the maximum stress is unaffected by the 




FIGURE 5. EFFECT OF FRICTION ON REVERSING SHEAR STRESS FOR RC TESTER 
(K = .8, v = .285, B-Model) 
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friction as is the depth to maximum stress (z/b - 0.38). There are some dif- 
ferences in the stresses near the surface, and it is impossible to be certain 
how these differences affect life. However, in the Lundberg-Palmgren theory, 
only the orthogonal shear reversal and its depth are important. If this 
theory is valid, it says that a traction coefficient of up to 0.1 should not 
alter life. There is, however, a difference in off-set between axial and 
tangential-friction, but this effect is not a part of fatigue theory. Fatigue 
tests are badly needed to clarify the effect of stress on life. 

The general observation from this aspect of the analyses is that the 
direction of the traction is not significant with regards to fatigue type stress 
for low traction coefficients typical of lubricated contacts assuming Lundberg- 
Palmgren's theory. It should be noted that at least three stress related 
factors that have not been included in these stress analyses could have a 
profound effect on fatigue. These factors are: 

(1) the effect of temperature (due to friction) on stress 

(2) the effect of non-Hertzian pressure on stress 

(3) the effect of asymmetric surface traction. 

All of these can produce an asymmetry in stresses which can amplify the stress 
reversal. In fact, there will be a different stress state in the driver and 
driven for effects (1) and (3). The influence of the non-Hertzian pressure 
should be the same on both bodies provided that we do not invoke hysteresis 
effects in the material. 

If we examine the influence of (1) and (3) on, say, the reversing 
orthogonal stress and assume that the traction is caused by the shearing of an 
elastic/plas tic-like material, then the surface traction in the latter half of 
the contact is larger than in the front. Also the traction stresses will be 
negative in direction on the driver and positive on the driven. This means 
that the stress reversal in the driven is Increased from the non-traction 
condition while it is decreased in the driver, possibly explaining some of the 
experimental findings reported elsewhere. The depth at which the maximum 
orthogonal acts is also reduced by an asymmetric traction stress. Traction 
produced by an elastic/plastic-like material may also explain why gears fall 
by pitting fatigue in the region directly below the pitch circle because the 
positive sliding conditions may produce a region of maximum surface traction 
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asymmetry and, hence, the largest magnitude of the reversing orthogonal shear 
stress. The thermal stresses produced by the sliding would also alter the 
stress profile for the orthogonal shear stress in a similar fashion but now 
because the driver (i.e., the faster moving body) would see a lower flash tem- 
perature than the driven. Thermal stresses would tend to increase the ortho- 
gonal shear stress in b th bodies but more in the driven than the driver. Also 
besides the increases in the amplitude of the stress reversal, the depth at 
which this occurs is reduced, giving a further reduction in fatigue life on 
the driver. Further research is needed in this area to establish the influ- 
ence of the effects just mentioned. 


DISCUSSION OF EFFECTS OF FRICTION ON FATIGUE STRESSES 


The subsurface stress model represents a useful tool for evaluating 
the role of friction on fatigue type stresses* The results of the evaluations 
of the stress analyses from the RC tester showed some interesting (although 
perhaps not profound) implications of the effect of lubricated traction on 
stresses* The purpose of this section is to further explore the effect of 
friction on stress* 

Figures 6 through 8 show predictions of the reversing orthogonal 
shear stress for various conditions of friction and aspect ratio (b/a). This 
reversing stress is normally considered the significant stress for 
fatigue* Note, however, in Figure 6, that friction does not alter the 
magnitude of the difference between maximum and minimum shear stress (i*e*, 
stress reversals). This implies that the magnitude of the orthogonal stress 
reversals are not dependent on friction* Changing aspect ratio does alter the 
magnitude of the stress as shown in Figures 9 and 10, but this dependence of 
stress on aspect ratio is not affected by friction* 

The effects of friction on the magnitude and depth of the maximum 
shear stress, are shown in Figures 11 and 12 for various values of aspect 
ratio* Some Increase in stress accompanies an Increase in friction. However, 
for fj ■ .1 and v » ,285 which are reasonable values for lubricated traction 
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FIGURE 7. EFFECT OF FRICTION ON REVERSING ORTHOGONAL SHEAR STRESS 
(K « .5, v « .285, z/b = .435 T-Model) 
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drives, only very small changes in stress depth and magnitude occur. Higher 
value of friction (f = .2) produces sizable increases in stress. The observa- 
tions concerning shear stress also occur with regard to the octahedral shear 
stress as shown in Figures 13 and 14. 

Figure 15 shows the effect of friction coefficient, f, on surface- 
shear-stresses-reversals for line (K = 0) contact. As can be observed, the 
magnitude of reversing stresses on the surface does not exceed the reversing 
orthogonal shear stress for values of f < 0.25. This implies that for this 
contact situation, fatigue-related stresses are not altered by friction for 
f < 0.25. Smith and Liu^) showed sizable effects at f = 0.33, which would 
be consistent with our computations. However, values of f of this magnitude 
should occur only in very poor lubrication situations and are not relevant to 
lubricated traction drive theory. As can be observed, friction does mate- 
rially alter the shear stresses near the surfaces, even for low (< 0.1) fric- 
tion coefficients. This implies that friction alters the overall stress field 
and could alter the propensity of the system to incur a surface or near sur- 
face initiated fatigue failure. Fatigue tests in conjunction with the stress 
model are badly needed to further extend our knowledge of the role of low 
levels of friction on fatigue. 

Figure 16 shows a shear stress plot of r xy versus t zx for the 
plane of maximum shear stress at the center of contact very near the surface. 
The magnitude of the stress reversals envelope (~ .6) is very consistent with 
the Smith-Liu^) computations and further illustrates the effect of severe 
friction (f » .33) on fatigue inducing stresses. However, the effect of 
friction on these stresses is not a factor at lower friction. 

In conclusion, it can be said that based on the computer models dis- 
cussed herein, low friction coefficient (f < .1) do not alter the fatigue 
stresses. At higher friction levels, changes will occur, but these cannot 
explain the observed effect of friction on fatigue. Other factors such as 
non-Hertzian pressures, thermal stress, and nonsymmetric tractions should then 
be Included to understand the roll of friction on fatigue stress. 







FIGURE 14. EFFECT OF FRICTION ON DEPTH OF MAXIMUM OCTAHEDRAL STRESS 
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FIGURE 15. EFFECT OF FRICTION ON MAXIMUM STRESS REVERSALS 
(K = 0, f A = 0, v = .285, B-Model) 
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SUMMARY 


The performance life of traction drives can be a very important 
factor in evaluating the overall effectiveness of the drive concept. Intui- 
tively, this life should be heavily related to the state-of-stress in the 
elements as a result of high loads and tractions. The purpose of the project 
has been to develop a mathematical tool for evaluating this stress-state. Two 
independent, but similar, approaches (one at Battelle and one by Tevaarwerk) 
to developing the stress model have been pursued. Essentially, the models 
involve the computations of subsurface shear stresses (and stress reversals) 
as a result of normal and frictional loadings. Where possible the models have 
been checked against existing theory and against each other. 

One use of the stress models has been to evaluate fatigue test 
methodology for a rolling contact (RC) fatigue tester. The standard RC tester 
consists of a small driven cylinder loaded between two larger crowned cylin- 
ders. The RC tester is a candidate device to evaluate the effect of traction 
on fatigue, either driving the crowned rollers (tangential slip) or by 
simply skewing the rollers (axial slip); the skewing technique would be 
easier. The stress model was used to evalute the effect of lubricated type 
tractions on stresses for a typical RC configuration. 

The general conclusion from the RC evaluations was that the fatigue 
inducing stresses (shear egress) reversals and depth of shear stress are not 
seriously altered by tractions (axial or tangential) which are typical of 
lubricated contacts. There are some changes in the stress field, but these 
are difficult to relate to fatigue. To further explore the effect of friction 
on stress, additional computer cases have been evaluated. 

The general conclusions from additional computer analyses is that 
friction coefficients of the order of 0.25 are required to materially alter 
the state of fatigue inducing stresses; this is much larger than would occur 
in lubricated contact. This implies then, that other effects not incorporated 
in the current stress-models should be included. These effects include non- 
Hertzian stresses distributions and thermal stresses. Both of these effects 
can impose nonsymmetrlcal (about the center of contact) shear stresses and 
thus will produce shear stress reversals which will superimpose on the normal 
and friction stresses. 
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APPENDIX A 


COMPUTATIONS OF THE STRESS TENSOR 


Stress Due to Point Loads 


The state of stress at any point on or beneath the surface of a 
system under concentrated contact is quite complex. To define this stress 
state at each point requires the computations of nine components of the stress 
tensor as related to contact pressures and tractions. The effect of pressures 
and tractions on stress were developed by Mindlin^^ for point loadings on the 
surface of a semiinfinite body using the coordinate system of Figure A-l, and 
lotting R^ =■ + y^ + z^. These equations appear as follows. 


Normal Load . 
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Hertzian Pressures and Surface Tractions 


For a point on the surface, the normal load could be written 
P N » pdxdy where p is the local pressure (A-19) 

Also, Pj * pf jdxdy and (.t-20) 

P A - pf A dxdy (A-21) 

where f^. and f A are the effective coefficients of friction acting trngential 
or axially to the surface. 

If the surface were subjected to a Hertzian pressure distribution, 
the pressure could be expressed (see Figure A-2): 


P 




(y - y j ) ? 

2 


(A- 2 2 ) 


where p^ is the maximum Hertz pressure, a and b are the major and minor *>xes 
of the contact ellipse, x and y are the coordinates with origins at the center 
of the ellipse and Xj and yj are locations of the point load relative to the 
center of the Hertz ellipse. The Battelle analyses were restricted to the 
case where - 0. Tevaarverk allowed y^ to be a variable and found y t - 0 to 
be the interesting region. 

The general approach for determining the nine components «/f the 
stress tensor at any point can be written 


'u 



dxdy , ( A-23) 


where represents nine coaponents of the stress tensor* for Herttian 
pressures end various frictional forces and are the results of Mlndlin's 
analyses (discussed in the previous section). In further coaputations, we 


. . .033. where Oj j is o x . 


* In tensor notation o^j iaplles Oj j, <Jj2» 0^3. o.j 
Oj2 in Tjjy. *22 ®y* ®23 *ys* 



A- 5 


shall use cr = a, /p such that all computations will be normalized on the 

ij ij H 

Hertz maximum pressures. 

It would serve no purpose to show all substitutions of Mindlin's 
equations into the general stress tensor equation. It can be noted that any 
of Mindlin’s equations containing the term y n (n odd) will produce a zero 
integral since the integration involves that function timer the even Hertz 
pressure over the whole pressure region for y^ =* 0. As an example of how the 


substitution occurs, consider the normal stress ° z for fj and f^ 


0 . 


a = 
z 
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where the lower limit L = - b Wl - ^ , 
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the upper limit U * + b y 1 ~ 2 
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FIGURE A- 2 


LOADING CHARACTERIZED BY HERTZIAN ELLIPSE 
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Scaling the Stress Tensor Equations 
The following substitutions were used to facilitate computations: 
z = z/b £ = x/z r\ = y/z K = b/a and x ^ = x^/b 
With these substitutions Equation A-24 appears 
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(A-25) 


where 


L = Xj/T - 1/z ij 1 - (zKn) 2 

U = Xj/Y + l/Y \ 1 - ( zKn) 2 , and 

— 2 2 

+ n +1 ♦ 


This type of substitution permitted reasonable computations near z * 0. 
In the computer solution 
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(otherwise) 
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All other stress componc were treated in the same manner using the 
general expression 




VT d£ d n 
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( A-27 ) 


where by substituting £ = x, n = y and z = 1 , all equations for o' become the 
^ ij 

equations for o’, . 


Review 


This is the most involved section in the report and may be difficult 

for the reader to readily grasp. Basically, all it involves is to take 

Mindlin's equations (A-l - A-18) , replace x with £, y with n, and z with 1 and 

integrate equation A-27. For example, suppose we want to calculate o' - » t 

l J x z 

from the normal load equation 


From the tangential 


From the axial load 
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(Equation A-5) 


(Equation A-ll) 


(Equation A-16) 


(A-28) 


(A-29) 


(A-30) 


R - 


where 


^1 + C 2 + n 2 



Equation A-27 then becomes 


TT 13 - -2 f f [c + £ 2 + Cn] Va dc dn . (A-31) 

“0 -=*> 

This equation can be solved by numerical computation. In the Battelle 
computations, the integral on £ was performed by a Gaussian quadrature 
approach and the y integration was performed by a Trapezoidal rule. 

Check of Numerical Solution 


Hamilton and Goodman ^ ^ developed an exact solution for circular 
contacts (K « 1). The numerical solution is compared with the exact solution 
in Table A— 1 ; the exact solution was generated as a part of the work of 
Tevaarwerk. The agreement tends to verify the numerical accuracy. 

Smith & Liu^ developed an exact solution for line contact (K * 

0). The numerical solution is compared with the exact solution in Table A-2. 


TABLE A-l. COMPARISON OF NUMERICAL SOLUTION WITH EXACT SOLUTION FOR K 
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APPENDIX B 


SOURCE LISTING FOR 
BATTELLE STRESS MODEL 






B-2 


21 

CONTINUE 

000730 


WRITEC1.1J) 

000740 

15 

FORMAT! / / . 20 X, * PLANE OF MAX SHEAR STRESS*/ / ) 

000760 


WRITE 1 1> 16) TMXXM.CNH 1, IXRFH).CN1(2. I XR FHI .CN1 ( 3. 1 XRFM ) » X 1 1 1 XRFM 

1 000760 


1>E13.5.* X-REF**t E13.5) 

000790 


WRITE (1.19) TAURVM, ASM (1,1 XRFM). ASM (2,1 XRFM), ASM (3. I XRFM) 

OOOSOO 

19 

FORMAT! * TAU RE VERS ING«* , E 1 3 . 5 » * AS < 1 ) • * . E 1 3. 5. • AS 1 2 I E 13. 5 , 

000810 


1 ♦ AS ( 3 ) # E 1 3 • 5 ) 

0C0820 


WRITE (It 16) 

000B30 

16 

F0RMAT(//»20X»* PLANE OF »AX OCT SHEAR STRESS ♦ //) 

000840 


WRITE 1 1.191 TMXXO»CN2(l»tXRFO),CN2(2,IXRFOI»CN2(3,lXRFQ),Xl(IXRFO 

t 000850 


11 ) 

000660 


WRITE (1.191 TAU* VO, AS (1. I XRFO ) . AS ( 2 . I XR F 0 ) , AS ( 3, I XRFO 1 

000870 


CALL RVSTRS(STRSS»CN3»A5»TRVM,TMAX,1,SIGN) 

000880 


WR ITE ( 1 1 17 ) 

000890 

17 

FORMAT (// » 20X. * RLANE OF REVERSING SHEAR STRESS *//» 

000900 


WRITE (1,18) TMAX»CN3(1,1 ) .CN3I2.1) , CN 3( 3,1) 

000910 


WRITE (1, 19) TRVM. AS(l.l) ,AS(2,1 ),AS(3, 1) 

000920 


WP ITE < 1 • 20 > 

000930 

20 

FORMAT ( 1H1 ) 

000940 


GO TO 8 

000950 

_ 7 

CONTINUE _ _ 



ENO 

000970 


SUBROUTINE STRSSC(IFLAG) 

000960 


COMMON/F/ IS.X1(22,1),GNU»Z»Y,XK,FT..A,MX1 

000990 


COMMON/C/ STRSS(3,3,22) 

001000 


DIMENSION FI ( 22t l 7 ).GI NT (22) » IC ALC < 22 > t TEMPI < 22 ) , FX< 22 ) t F IY( 1 7) 

001010 


1 « Y 1 ( 17 ) 



PIO « 141 59 

001030 


MY-17 

001040 


M X »M X 1 

001050 


CALL XPOS 

001060 


TFIIFLAG.EO.lt GO TO 10 

001070 


00 9 ISTRSS-1.6 



1S-ISTRS5 

001090 


00 4 JY-ltMY 

0011C0 


IF(JY.EO.l) Y»0. 

001110 


IFIJY.E0.2) Y-.2 

001120 


IF ( J Y.EO. 3 ) Y-.4 

001130 


IFUY.EQ.4) Y-.6 _ 

Fnnrw 


IF (JY.EQ.5 ) Y*.B 

0011)0 


IF ( J Y . 6 0* 6 ) Y*1 • 2 

001160 


IF < JY.EQ.7I Y • 1 . 6 

001170 


IFUY.FQ.8> Y«2 • 

C011 60 


IF ( J Y.E 3*9 ) Y-2.5 

001190 



FTrrmHW 


IFUY.EQ.11I Y*4* 

001210 


IF(JY*E0,12) Y*6. 

001220 


IF(JY»E0*13) Y*8« 

001230 


IFUY.E0.14) Y-10. 

0012*0 


IF UY. EO. 15) Y*20 • 

001250 





IF ( JY.F0.17) Y-50. 

001270 


Y1UY)-Y ....... . . . . - 

001200 


61 — z.n 

001290 


IHBl.lT.-20.) BL— 20. 

001)00 


BU— 61 

001310 


1ELE»1 



CALL GAUSS(NX»l*20»l»l.E-3.1.E-3.BL.BU.GINT.lCALC»TEm,FXI 

001330 


DO 5 IX*1» MX 

001340 


FI(IX,JT)»GINT(lX) 

001350 

5 

CONTINUE 

001360 

4 

CONTINUE 

001370 


B-3 



on 6 i*>i,mx 

001380 


on 8 ir*i,"T 

001390 

6 

Fircm*FiMx»m 

001400 


£AkL _Y I NTtfir.fTl»SlGliHT) 

001410 


IFIIS.EO.il STRSSll.l. IXI-SIG1 

0)1420 


IFdS.EQ.2) $TRSS(2,2,IX)-SIG1 

001430 


IF€IS.E0,3I STRSS (3# 3» IX )■$ IG1 

001440 


IF* IS.E3. 4 1 STRSSU»2»IX|.$TBS$C2»l,m*SlGl 

001450 


IFdS.EQ.5l STRSS (2, 3. I X I -STR SS ( 3, 2 . 1 X 1 • S I G 1 

001460 


IHIS.E0.6) ST«SSU»3»I»fSTRSSO»l«lKl»S161 

001470 

6 

CONTINUE 

001480 

9 

CONTINUE 

001490 


RETURN 

001500 

10 

CONTINUE 

001510 


DO 12 ISTRSS-1,6 

001520 


1S«ISTPSS 

001530 


00 11 I * ■ 1 » N X 

001540 


X»Xl(IX#i> 

C01550 


XK1»I1 .♦! » **2»Z**2 

001560 


XK2-d.-X)4R2d*R2 

001)70 


C1*SQRT(XK2/XK1) 

C01580 


CZ»S0RT(2.*C1HX<1*xk2-4.)/ik1) 



SI*PI/XK1*(1.-C1I/CI/C2 

001600 


SI8-PI/XK1** l.*Cl 3/C1/C2 

001610 


SXN»-l,/PI*Z*dl.tZ.*X**Z»2.*Z**2)*SlB-2.*Rl-3.*X*SI I 

001620 


SZN»-Z/PIMSIB-X*SII 

001630 


T*ZN— l./PI#Z*R2RSI 

001640 


nio(ruuiiw«mw«u}4uwifljuj ufiiuuw nm 



TXZT-FT*CTZ 

001660 


TYZT-FA*CTZ 

001670 


CXT— l./PI** ( 2. »X **2-2. -3. *Z**2)* SI *2. *Pl*X *2. *<1.-X**2-Z**2)*X 

001660 


1 ♦ S I B 1 

001690 


SXfFMCXT 

001700 


$tt«fa*cxt 



SZT»-<FT*FM/PI*Z**2«SI 

001720 


IF(IS.EQ«1I STRSS U»l, IXI.SXNPSXT 

001730 


IF ( IS.EQ.2 ) STRSS (2.2* IX)»GNU*($XN.SXT»SZN»SZTI*5YT 

001740 


IF(IS.E0.3) STRSS<3,3. IXI-SZNMZT 

001750 


IF ( IS. E0.4 I STRSS d, 2, IX I -STRSS (2,1, IX) >0. 

001760 





IFCIS.E0.fr) STRSS(1,3,IX)«STRSS(3,1 »IX)-TXZNpTXZT 

001780 

11 

CONTINUE 

001790 

13 

CONTINUE 

001800 


RETURN 

001810 


ENO 

001820 





DIMENSION FIC17I.Y1U7) 

001840 


C1«FI(1I«( YH2l-Tlim 

001850 


C2»FUHn*iniNY»-VHNY-in 

001860 


SUN.O. 

001870 


MYl-NY-1 

001880 



001890 

1 

suN.suN»Fim*(Yid»n-rid-in 

001900 



SIG1»C1*C2*SUN . . _ 

— 001910 


RETURN 

001920 


END 

001930 


SUBROUTINE FEVM.IX.FI, JF.ICUCI 

001940 





DIMENSION FX(22)» ICALC (22) 

001960 


RIO. 141)9 

001970 


X2>X*«2 

001990 


rjiYMZ 

001990 


GNUl«l.-2.RGNU 

002000 


B-4 



R»$amx**2*Y**2*i . ) 

R2»Q*#2 

R3*R^3 

R5«R**5 

002010 

002020 

002030 

002040 


RZ-R*1. 

002030 


IF(ISTRSS.EO.l) 1,2 

002060 

1 

Cl«<GNUl/*3-3'M2/R5-GNUl/R/ftZ*(l'-X2/RrRZ-X2/*2l > / ? • / P I 

002070 


C2 • (GNU1/R 1- 3. RX2 /R 3-CNU1/R/ R Z**2* 1 3.-Y2* ( 3.*R»l.l/R2/RZll/2./PI**C02060 


C3«0. 

C02090 


GO TO 10 

002100 

2 

I F t ISTRSS «£0 .2) 3»* 

002110 

3 

C1«<GNUWR3-3.*Y2/R5-GNU1/R/RZM1.-Y2ZR/RZ-Y2/R2» 1/2./RI 

002120 


C2«IGNU1/R 3-3. ♦Y2/R5-GNU1/R/R 2**2* U.-Y2*t3.*R*l. I /R2/ R Z ) 1 *X / 2. 

/ P 10021 30 


C3-0. 

002150 


GO TO 10 

002150 

4 

IFUSTRSS.EQ.3) 5»6 

002160 

5 

Cl*-3./R5/2. /PI 

002170 


C2*-3*/2./Pl*X/R5 

002180 


C 3*0 • 

002140 


GO TO 10 

002200 

6 

IFUSTPSS.EQ.4) 7,8 

002210 

7 

Ci*0. __ 

. 002220 


C2-0. 

002230 


C3»<-3*4<2/R5-GNUl/«/RZ442Ml*-X2*(3,PR4l.)/R2/RZn/Z./PlM 

002250 


GO TO 10 

002250 

e 

1 F ( ISTRSS « EQ • 5 ) 9,11 

C02260 

9 

C1«0. 

002270 


C 2 »0 • 

_ClO22B0 __ 


C3 — 3./2.JPMY2/R3 

002290 


GO TO 10 

002300 

u 

Cl»-3*/2*/PlM/R3 

002310 


C?»-3*/2*/PI*X2/R5 

002320 


C3-0* 

002330 

10 

AR G»C1*FT_*C2 ♦FJ^C 3 

002350 _ 


00 15 I • 1 , JF 

002330 


IFUCALC(mi5,13 # 13 

002360 

13 

CONTINUE 

002370 


ARGU-l'-<2B*X-U<Itin**2-(2B*Y*XK)*t2 

0023B0 


IF ( ARGU* L T .0 • 1 13,16 

002390 

15 

FYm.i.i-i2_ 

002400 


60 TO 17 

002410 

15 

CONTINUE 

002420 


f I<l>«SORTURCUI***G 

002430 

17 

CONTINUE 

002440 

15 

CONTINUE 

002430 


RETURN _ ... _ 

- 002460 


ENO 

002470 


SUBROUTINE M INC ISTRSS# IX* SH» SIN# SHNI 

002480 


OINENStON STRSSI 3.3*22 1 

002490 


Cl"STR$SU.l.ll>»STRSS<2*2»lXt«STR$St3»3»m 

002300 


C1«-C1 

002310 


C2»STRSSC1»1* IR1.RSTRSS I2»2*J*J»STRSS11»1» IRLPSTRSSJl* 1» IX1> _ 

001320 


lSTRSS(2.2.in*STRSS(3.3.m-ST«SS(1.2,H»*«2- 

002310 


2STRSSU.J. I*)**2-STRSS(2j3.UI**2 

0023*0 


C!«STRSStl#l#IM>*ST»SS<2.2,in*ST»SS<3.3.m* 

002330 


12.*STRSS(1.2tin*STRSS({j)»in*STRSS(l»3f 1*1- ... 

0023*0 


2ST»SS< 1.1# tn*ST*SS <2.3.1* )*«2- 

002370 


3STMSSi2.2.m*STRSS(l.l.m**2- _ . _ 

0023B0 


4$TISS<3,3, m*STRSSU»2,IX>**2 

002390 


C3»-C1 

002600 


CAU CUBIC IC1,C2,CJ, SI, SZrSJ) 

002610 


S«»Si 

002620 


IF(S2'6T.$*I SN«32 

002630 



B-5 



IMS3.GT.SH) SM-S3 

002640 


snn.si 

002630 


IMS2.lT.SnN) SNN-S2 

002660 


:F(S3.1T,SMN) $NN«$3 

002670 


Sin-Si 

0026B0 


IF (S2.GT.S HN.ANO. S2.LT. SM) SIN. $2 

002690 


IF (S3.GT.SnN.ANO.S3.LT.Sn) Sin. S3 

002700 


RETURN 

002710 


ENO 

002720 


SUBROUTINE CU3ICICl»C2fC3#Sl*S2#S3) 

002730 


pio.ms 1 ) 

002760 


*«1./3.*I3.*C2-C1*R2) 

002750 


B.1./27.«(2.*CI**3-9.*C1*C2»27.«C3I 

002760 


IF( A.LT.O. ) GO TO 3 

002770 


S1-S2.S3-0. 

002780 


RETURN 

002790 

3 

CONTINUE 

002600 


PHI.*C0S(-B/2./$3RT(-4**3/27. ) ) 

002910 


$l«2.*SQRT(-A/3. l*CQ$(PHI/3. t-Cl/3. 

002920 


S2.2.*S3RT(-*/3, ) *CQS I RHI / 3 . »2 .*P I / 3 . I-Cl/J. 

002930 


S3«2.*S0RT(-*/3.)PCDSIPHI/3.»N.PPI/3. )-Cl/3. 

002860 


RETURN 



ENO 
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SUBROUTINE *N51PI$TRSS»S1»S2»S3*BN#IK»S16NA) 

002870 
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002880 


SIl-STBSSI 1,1*IX) 

002800 
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002900 
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002910 
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002920 
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002930 
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002960 
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002930 
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002990 
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003000 

3 
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003020 


GO TO 1* _ _ 
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003060 
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003090 
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003060 


CO TO 14 

003070 

6 
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7 

BNtl.lfl. . 



RNII,1)«*NII,2).0. 

003100 


CO TO 1* 

003110 

• 

IM*BSfS12).lT.E*.*N9.tBSm3>.lT.CR) 9.10 

003120 

1C 

A— 1512*513-1 SU-SI*S23)mi2*»2-<Sn-S)*<S22-S>) 

001130 


CliMI-S12*A-Ull/tSll-SI)**l*l.*A9*2 

003160 





AN(I»2)«A*ANtI»3> 

003160 


BN 1 U 1 1 • ( -512*AN LI* l l-UJ* ANUjll 1/lUl-l L_ 

.003170 


60 TO 16 

001100 

9 

IHAISt$l2 )tlT.i*«ANO. AlSfS23>tlT • !*) 11#12 

_ _ 003190 

n 

A— CS12*S23-CS22-S)*S13)/(S23**2-CS22-S)ECSJJ-S) > 

003200 



003210 


ANII»1)>SI0NA*S0BT(1./C11I 

009820 


ANII.3) S A*ANII»1) 

003230 


ANU.2).(-S2)*ANI J.3)-S12*ANU.1))/(S22-S> 

003260 


CO TO IB 

009290 

12 

CONTINUE 
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A.-<S13*S23-S12*< S33-S M/IS1 3**2- < S ll-S > • < $ 33-S > > 


003270 



CU-I <-S23-S13*A)/IS33-S ) )«*2*1 .♦*♦*2 


OC 32 80 



ANII.2)«SIGNA*S3RT(1./C11> 


003290 



AN(I,l)-A*Afi(I»2) 


003300 



ANUt3>«<-Sl3«AN| 1 , 1) -$ 2 3* AN ( 1 , 2 ) )/( $ 3 3-S ) 


003310 


14 

C ON T I NUF 


003320 


2 

CONTINUE 
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003340 
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SUBROUTINE * VST PS (STPSS*AN#AS# TP VN* T"U X « I » S I GN ) 
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COHHON/F/ ISTRSS. *1122.1 ) . GNU. 78. T . << .F T . F A , AX 


003370 



01 AENSION ST R SSI 3,3.22). AN < 3.221. A$( 3,22) 


003380 



trva-o. 


003390 



TAXX«0. 


003400 



ER»l.E-3 


003410 



00 1 I*S-1 »21 


<1014 20 



AS1»-1.*.1*F10ATI IAS-1 ) 


003430 



IF<A8S<AN< 2, I M.IT.ER.AN0.A3$(AN(1> I > ).IT.ER) 8,9 


003440 


8 

AS3-0. 


003430 



AS2«SIGN«S0RT(1.-AS1**2» 


003460 



GO TO 13 


003470 


9 

.F(AnS(AN(l>I)).LT.ER.AN0.ABS(AN(3#I)).LT.ER) lOill 


__M148Q . 


10 

AS 2 *0 • 


003490 



AS3«SIGN*S09T(1. -»ASl*t2) 


003300 



GO TO 13 


003310 


n 

IFCA65(An(2# I) ).IT.EP.AND.AB$< ANH# IIKIT.EP) GO TC 12 


003320 



IF(A8S(AN(2fl))*lT*£P) 16tl7 


003330 


16 

AS3"-ANtl* ntASl/ANOnI) 


003340 



CC-1.-AS1**2-AS3*R2 


003330 



IFICC.LT. 0.1 63 TO 22 


003360 



AS2«SIGN»50RT (CC 1 


003370 



GG TO 13 


003380 


17 

t c <A8S<AN<3,I)).LT.ER> 18.19 
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16 

AS2 — AS1*AN< 1. I)/AN<2, I ) 
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CC*1.-AS1A#2-aS2AA2 
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IF ( CC . L T « 0 • ) GO TO 22 


003620 



AS3*S ICNASOP T ICC > 


003630 



GO TO 13 


003640 


19 

tHAS$UN(l»n)*lT.m 20*21 


003630 
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as3»-an<2, n/AN(3.n*A$2 
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GO TO 13 
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21 

A>AN<2.I)«*2«AN<1,1)«*2 


003690 
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003700 
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003710 
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IF ICC .LT.O . ) 60 TO 22 


003730 
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AS2«-I A$3* AN< 3. I 1 AASl* ANI1.II1/4NI2.11 


003730 


13 

COUTINUI 
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001770 
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001790 
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_ 

001900 



l«AN<).!)*iS3*STRSS<).).m*ANfl.l)*4S2*SmS<l,2.!l> 


0011 10 



2»A*l2,IJRASl*$TRJSU,2/m*AMtl,n*»$J*JT«JStl, 1» in 


001920 



]*tNUin*tSlASTItJ(lt)iIll»ANU.n*A3)«STIJSIMiIII 


001910 








IPtTAU.GT.TNkl 3# 4 


001990 


) 

TNI»TAU 


001960 



AST1.AS1 


001970 



AST2*AS2 


001990 



AST3*AS3 


001990 
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IMT4U.lT.THN) 9.6 


003900 

5 

TNN»T4U 


003910 

6 

CONTINUE 


001920 

.2 

CONTINUE _ 


003930 


TRV-THI-THN 


003940 


IF( ABS<T*tf I.GT.ASSmVH) ) TRVM-TRV 


00 3 Q 50 
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003970 

23 
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22 

CONTINUE 
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CONTINUE 
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RETURN 
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12 
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15 

T0BH4TC • SOLUTION NOT A RRROBR I AT E ♦) 


OOA060 


BE TURN 
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END 
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SUBROUTINE ANGTRAN(AN.BN'CNiIX) 


C0A090 


OIHENSION AN(3«3l»8Nm»CN(3»22) 


GOAIOO 


CN(1»IX)*5N(1)AAN(1.1).BN(2)RAN(2.1).9N(3)RAN(3»11 
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CNI2*IX)»8MI)AAN<1#2)»6N(2)4AN(2*2)*8N(3)RANI3*2) 


00A120 
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RETURN 
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SUBROUTINE XROS 


00 A 160 
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OOA190 
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00A20O 
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00A210 
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0OA22O 
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004240 
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004230 


RETURN 


0OA260 

2 

IFMX.EO.ll) 3. A 


O0A27O 

3 
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00A290 
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CO A 300 
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00A31O 
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0OA12O 
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00A31O 
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OOA3AO 




004330 
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0OA36O 


*1(10, 1M. 9 


00A370 
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RETURN 


00A3R0 
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5 
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_ 

006630 
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■ miii 
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006680 
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006910 
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RETURN 

004530 

004540 

004550 
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PRINT 10 
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10 

FORMAT! * MATRIX NOT AVAILABLE * ) 
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RETURN 
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END 
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